Quantum transport modeling of Fe/MgO/Fe magnetic tunnel junction with FeOo.5 

buffer layer: the effects of correlations 



(N 

o 

(N 
O 

Q 



CZ2 



o 
o 



> 

(N 



X 



V. TimoshevskiiQ Yibin Hu, E. Marcotte, and Hong Guo 
Department of Physics, McGill University, 3600 rue University, Montreal, Quebec, Canada H3A 2T8 

(Dated: December 7, 2012) 

We report ab initio simulations of quantum transport properties of Fe/MgO/Fe trilayer structures 
with FeOo.5 buffer iron oxide layer, where on-site Coulomb interaction is explicitly taken into account 
by local density approximation + Hubbard U approach. We show that on-site Coulomb repulsion in 
the iron-oxygen layer can cause a dramatic drop of the tunnel magnetoresistance of the system. We 
present an understanding of microscopic details of this phenomenon, connecting it to localization 
of the Fermi electrons of particular symmetry, which takes place in the buffer Fe-O layer, when 
on-site Coulomb repulsion is introduced. We further study the possible influence of the symmetry 
reduction in the buffer Fe-0 layer on the transport properties of the Fe/MgO/Fe interface. 

PACS numbers: 85.35.-p, 72.25.-b, 85.65.+h, 73.20.-r 



I. INTRODUCTION 

Due to its high tunnel magnetoresistance (TMR), 
studying of Fe/MgO/Fe magnetic tunnel junctions 
(MTJ) has evolved into an active research area dur- 
ing last decade. Indeed, the TMR value of 180- 
250%, reported by several research groups for MgO-based 
systems^, is the highest TMR, experimentally measured 
at room temperature. These high values of tunnel mag- 
netoresistance make the MgO-based MTJs excellent can- 
didates for spintronics applications: magnetic random ac- 
cess memory^, programmable logic elements^, and mag- 
netic read sensors. 

Since the original theoretical prediction 5 of large TMR 
effect in Fe/MgO/Fe interfaces, atomic-level simulations 
played a key role in understanding of unusually high 
TMR values in these systems, and a lot of progress 
has been achieved in this field. Butler et ai 5 -, using a 
Korringa-Kohn-Rostoker (KKR) technique, has provided 
an elegant explanation of physics behind the Fe/MgO/Fe 
magnetoresistance, connecting it with the match of the 
symmetry of the Bloch states in the electrodes with 
the symmetry of the evanescent states of the MgO bar- 
rier. Later, Waldron et al 6 , employing Non-Equilibrium 
Green's Function (NEGF) formalism, combined with 
Density Functional Theory (DFT), explained a voltage 




FIG. 1: Atomic structure of the two-probe Fe/MgO/Fe de- 
vice with Feo.sO buffer layer. The system is periodic in the 
directions (X, Y) parallel to the interface. The leads are semi- 
infinite in Z-direction. 



dependancc of TMR in Fe/MgO/Fe interfaces. However, 
there is still an important issue which is the quantitative 
understanding of physics behind Fe/MgO/Fe system: the 
absolute values of experimentally observed TMR ratios 
cannot be reproduced by ab initio simulations if assum- 
ing ideal device models: these ideal models predicts much 
larger TMR than that has been observed. There may be 
a number of reasons behind the theory /experimental dis- 
crepancy such as oxygen vacancies in the MgO layeri^, 
the existence of a FeO buffer monolayer^S, other atomic 
defects^, and interface roughness^. These structural im- 
perfections are clearly detrimental but may be overcome 
by improving the material quality and the fabrication 
process of the devices. It is the purpose of this work to 
investigate a possible physical reason - as opposed to the 
structural reason, for the reduced TMR value as observed 
experimentally. 

In particular, we theoretically investigate an impor- 
tant piece of physics that is related to the presence of 
a FeO buffer monolayer formed as a result of electrode 
oxygenation during the MgO deposition 9 . Calculations 
already demonstrated^ a reduction of TMR from that 
of the ideal device structure in the presence of the FeO 
layer, although even at 60% oxygen concentration in the 
layer, the predicted TMR is still much greater than that 
observed experimentally. On the other hand, it is well 
known that bulk FeO crystal is a Mott-Hubbard insula- 
tor whose electronic structure cannot be correctly pre- 
dicted by Local Spin Density Approximation (LSDA) 
commonly used in density functional theory (DFT) cal- 
culations. This is due to an incorrect treatment of strong 
correlation in the FeO compound, as such LSDA predicts 
FeO to be metallic while in reality it is an insulator with 
a well-developed band gap. It is therefore very inter- 
esting to ask if a strong correlation exists in the single 
monolayer of FeOi-a; that is at the interface of between 
Fe and MgO in the Fe/MgO/Fe trilayer so as to influ- 
ence TMR. If it does, a localization of the Fe 3d-states 
is expected that may lead to a substantial drop of trans- 
mission coefficients for both up- and down-spin channels 
and, as a result, to a significant change of TMR. To our 
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knowledge, the only attempt to take into account corre- 
lation effects in FeO layer in Fe/MgO/Fe interface was 
taken by Mirhosseini et aU&, who applied self-interaction 
correction (SIC) to LSDA to calculate spin conductance 
and magnetoresistance of Fe/FeO/MgO/FeO/Fe. These 
authors obtained a significant increase of the binding en- 
ergies of Fe 3d-states in the "spin-up" channel in FeO 
layer. As a result, it was found that SIC-corrected treat- 
ment moderately reduces the TMR ratio at energies be- 
low the Fermi level (Ep), from where the Fe 3c? majority 
states are removed, while a 10% increase of TMR at Ep 
and above is observed. Since correlation effects may be 
important and it is a very complicated problem, further 
investigation along this line is warranted. 

Here we present results of our first principles calcu- 
lations of quantum transport properties in Fe/MgO/Fe 
interface with FeOo.s iron-oxide separation layer with 
the oxygen content close to the experimental one (60%). 
We use a state-of-the-art quantum transport technique 
based on density functional theory (DFT) combined 
with Keldysh nonequilibrium Green's function (NEGF) 
formalism^ 3 -. The correlation effects in FeOo.s layer are 
taken into account within the LDA+[7 approach which 
explicitly adds the on-site Coulomb repulsion energy (U) 
for Fe 3d-states directly to the DFT Hamiltonian. The 
value of U is calculated from first principles for a partic- 
ular atomic structure of the interface using a constrained 
LDA technique^. Our results demonstrate that the con- 
sideration of correlation effects in the FeOi- x layer leads 
to a dramatic drop of magnetoresistance in Fe/MgO/Fe 
interface, and brings it to the values comparable to that 
observed in experiments. The microscopic details of this 
substantial TMR reduction can be understood by analyz- 
ing the effects of electron localization on scattering states 
of different symmetry. 



II. CALCULATION DETAILS 

A. Interface structure 

We consider a model Fe/MgO/Fe interface structure 
with five layers of magnesium oxide stacked between two 
semi-infinite Fe-bcc electrodes, oriented in the (100) di- 
rection (Fig. Q]). The system is periodic in the plane 
transverse to the interface direction. The unit cell vec- 
tors in this plane are rotated 45° with respect to Fe-bcc 
unit cell, and their lengths equal where a is the cell 
parameter of Fe-bcc lattice. This choice of the system 
geometry allowed us to model the interface FeO mono- 
layer on each side of MgO slab with 50% oxygen occu- 
pation (Fig. [T]), which is close to 60%, observed experi- 
mentally by Meyerheim et al-. The structure geometry 
was completely relaxed with respect to atomic positions 
and intcrlayer distances using highly accurate all-electron 
Linearized Augmented Plane Wave Method (LAPW) as 
implemented in Wien2k 15 program package. We paid 
special attention to geometrical properties of the inter- 



face, as the earlier study of Waldron et al demonstrated 
a substantial change in TMR, when the interface atoms 
are displaced by only 1% of their bond length 6 . Our the- 
oretical optimized geometry showed a good agreement 
with results of X-ray diffraction experiments performed 
for the MTJs having 60% oxygen content of the inter- 
face FeO layer—. In particular, the obtained Fe-0 dis- 
tances within the interface FeO layer and between FeO 
and MgO layers were found to be 2. 03 A and 2. 38 A re- 
spectively, which are in very good agreement with exper- 
imental data (2.03A and 2.35A). 

Having obtained the relaxed geometry of the interface, 
we proceed to calculate the electronic transport proper- 
ties. We used the MatDCal program package^ which is a 
DFT-based code, employing a state-of-the-art quantum 
transport technique based on the Keldysh nonequilib- 
rium Green's function (NEGF) formalism. The method 
uses standard norm-conserving pseudopotentials— and 
an s, p, d double-zeta LCAO-type basis set. The NEGF- 
DFT formalism is excellently suited for self-consistent 
electronic transport calculations of multi-atom inter- 
faces with non-periodic (open) boundary conditions. A 
detailed description of theoretical base of the method 
can be found in our previous publications 13 ' 18 . We 
found a 20 x 20 K-pomt mesh to be sufficient for self- 
consistent calculations of the two-probe device. Much 
denser (10 3 x 10 3 ) sampling was further used for calcu- 
lation of transmission coefficient T. A TMR ratio was 
calculated at zero bias voltage using the following for- 
mula: TMR = (T-|~|- — T-fi)/T^i, where Tff and 
correspond to transmission coefficients calculated at the 
Fermi energy for parallel and anti-parallel spin configu- 
rations of iron leads. 



B. Correlation effects 

To study the effects of electron localization on elec- 
tronic transport, we went beyond the standard NEGF- 
DFT formalism, and implemented the LDA+ U scheme 
into the MatDCal code. We followed the original works of 
Anisimov et a ? 14 ' 19 i 20 , and details of our implementation 
of the LDA+ U method within NEGF formalism is sum- 
marized in the Appendix. According to the philosophy of 
LDA+ U approach, the atomic system is divided into the 
weakly-correlated part where the exchange-correlation is 
treated within LDA; and a strongly-correlated part where 
the Coulomb (U) and exchange (J) energies of interac- 
tion of two electrons on the same atom are explicitly 
taken into account. In the case of oxidated Fe/MgO/Fe 
system, we treat the interface FeOo.5 layer where we ex- 
pect the electron localization to play a major role as a 
highly correlated subsystem and apply LDA + U treat- 
ment to the Fe-3d states of this layer. 

To correctly apply the LDA + U method, proper val- 
ues of U and J parameters need to be supplied. These 
parameters are usually either taken from experimental 
data, or a series of calculations are performed for dif- 
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ferent values of U and J to check how the results de- 
pend on them. In our case, we have a single monolayer 
of FeOo.5 sandwiched between an Fe electrode and the 
MgO barrier. To our knowledge, no experimental data 
about the on-site Coulomb repulsion energy in this type 
of interfaces have been reported, hence we have calcu- 
lated it from first principles. Note that the U parameter, 
defined as the cost in Coulomb energy to place two elec- 
trons at the same sitei^, can be calculated by following 
the original ideas of Anisimov and Gunnarsson 14 within 
their particular implementation of the LAPW method, 
proposed by Madsen and Novak^i. In our calculations, a 
supercell parallel to the interface plane was constructed, 
and the hybridization of 3ci-states of one of the Fe atoms 
in the FeOo.5 layer with all other orbitals was switched 
off. By performing a series of self-consistent calculations 
for this system with different values of occupation num- 
bers of this artificially atomized 3d-shell, we calculated 
the U e f f — U — J parameter as a difference between the 
corresponding eigenvalues^ 1 . This way, we found that the 
effective Coulomb repulsion for Fe 3d-electrons in the in- 
terface FeOo.5 layer is about 4.6 eV. This value is some- 
what smaller than the one calculated by Anisimov et al 
for bulk FeO crystal (5.9 eV)^. The reduction of U in 
our system is expected due to the presence of the Fe elec- 
trode on one side of the FeOo.5 layer: the free electrons 
in the Fe electrode increase the screening of the on-site 
Coulomb repulsion in the FeOo.5- 



of magnitude as the results obtained in other similar 
calculations^, as well as the LDA results obtained for 
oxigenated Fe/MgO/Fe interfaces (~ 1500%}^. As soon 
as Coulomb repulsion of 4.6 eV is taken into account, the 
TMR ratio drops dramatically down to 60%. To trace 
the origin of this drop, we analyze the transmission con- 
tributions to TMR from different spin channels for PC 
and APC (Tabled]). We notice that the transmission in 
APC for both spin channels reduces by a factor of ~ 3 
when Coulomb repulsion is introduced. Even smaller is 
a reduction in spin-down channel for both PC and APC 
(factor of ~ 1.7 — 2.7). A really defining drop in trans- 
mission is observed in the spin-up channel for the case of 
PC. In this case the transmission coefficient drops by a 
factor of 46, which constitutes a principal reason behind 
the dramatic reduction of TMR. 

We should note that our LDA+J7 results are some- 
what opposite to those of Mirhosseini et aU^, who used 
the SIC-corrected LDA approach. We clearly see a signif- 
icant drop of TMR ratio upon "turning on" correlations 
in FeOo.5 layer, while Mirhosseini et al did not observe 
significant changes of this quantity at energies equal or 
above the Fermi level. Further, our results demonstrate 
a principle reason of this drop - a dramatic reduction of 
transmission in the "spin-up" channel for PC, while the 
SIC correction mostly reduces the APC current. 



B. Scattering states 



III. RESULTS AND DISCUSSION 
A. TMR ratio 

First, we performed full self-consistent calculations of 
the interface for Parallel (PC) and Anti-Parallel (APC) 
spin configurations of Fe leads for the cases of f/ e // = 0, 
as well as with "turned on" Coulomb repulsion. Further, 
for the case of each U -value, a contribution to transmis- 
sion coefficient from spin-up and spin-down channels has 
been obtained, as well as a final value of TMi?-ratio. The 
obtained results are summarized in Table |H The most 
striking result is the dramatic drop of the TMR value 
for as soon as on-site Coulomb repulsion is taken into 
account for Fe 3<i-electrons in the FeOo.5 interface layer. 

For a pure LSDA calculation (U e ff — 0), we ob- 
tain a TMR of 1942%, which is of the same order 



TABLE I: Calculated contributions to transmission coefficient 
(xlO -6 ) from different spin channels for the cases of U — 
and U = 4.6 eV. The calculated TMR ratio is also provided. 



Configuration 


T t 


n 


TMR 


Ueff 


PC 
APC 


95.57 
2.56 


0.81 
2.16 


1942% 


0.0 eV 


PC 
APC 


2.08 
0.80 


0.47 
0.80 


60% 


4.6 eV 



To establish a more detailed physical picture, it is in- 
teresting to follow how the absolute square of the tun- 
neling electron wavefunction (scattering state) changes 
along the tunneling direction. As the drop in TMR ra- 
tio is caused by substantial reduction of transmission in 
spin-up channel for PC, we calculate the scattering state 
for this particular spin channel in PC. Only states at the 
T-point and at the Fermi energy are considered. 

In their pioneering theoretical study, Butler et a& 
demonstrated that tunneling in the spin-up channel al- 
most exclusively takes place through the state of Ai sym- 
metry, while other bands, crossing the Fermi level of the 
Fe lead, rapidly decay in the MgO barrier. In our system 
we also observe two low-decay scattering states with sim- 
ilar profile, while others decay faster in the MgO barrier. 
The presence of two of these states, as compared to But- 
ler's setup where there is only one state of Ai symmetry, 
is explained by the choice of our unit cell, where there 
are two Fe atoms per layer in the direction perpendicular 
to tunneling. We plot the low-decay scattering states in 
Fig. [2j It can be noticed that these two states are sub- 
stantially different. Even without the Coulomb repulsion 
U, there is a factor of 10 3 difference in the squared ampli- 
tude of these states after tunneling through the barrier. 

To get insight into the physical origin of this differ- 
ence, we calculate the spectral composition of the two 
scattering states, shown in figure [H by projecting them 
on atomic orbitals of the interface Fe atoms. We were 
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FIG. 2: Distribution of squared electronic wavefunction (scat- 
tering state) across the interface for two low-decay tunneling 
states of Fe/MgO/Fe interface with FeOo.5 buffer layer. The 
wavefunction is integrated in XY plane. Calculations are 
done at the Fermi energy, for spin-up channel and at the T- 
point only. The arrows mark the positions of the FeOo.5 layer 
on both sides of the interface. 



surprised to find that the orbitals, forming the Ai-state 
(s,d z 2) in the system without FeO layer (or with 100%- 
populated FeO layer), were not the only ones forming 
these low-decay states in the case of our structure with 
FeOo.5 iron oxide layer. Additionally, a significant contri- 
bution of d x 2_ y 2 Fe states was found in both scattering 
states, presented in Fig|2] In fact, this contribution is 
comparable in value to the one of d z 2 orbital, and even 
exceeds it for one of the scattering states. The calcu- 
lations with U = showed that the contribution of d z 2 
and d x 2 +y 2 orbitals is respectively equal to 54% and 30% 
for scattering state (a), and 19% and 56% for scattering 
state (b), FigO We see that the scattering state (a) is 
mainly formed by the d z 2 states, while the d x 2 +y 2 orbitals 
mostly contribute to the formation of scattering state (b) . 
This difference in electronic structure is the underlying 
reason of the significantly different decay rate of these 
two scattering states. A much lower tunneling amplitude 
of state (b) can be understood as follows. The d x 2 +y 2 
orbitals, which contribute mostly to this state, are ori- 
ented in the XY plane (parallel to the interface), which 
reduces their overlap with p z -orbitals of MgO layer. At 
the same time, the d z 2 orbitals, which dominate in the 
formation of scattering state (a), are oriented along the 
Z axis (perpendicular to the interface). This orienta- 
tion increases their overlap with MgO oxygen states and 
therefore creates conditions for higher tunneling current. 



C. Electron localization 

Let us now concentrate on the scattering state Fig[5^a) 
which has the lowest decay rate, and therefore mostly 
contributes to the tunneling process. We see that "turn- 
ing on" Coulomb repulsion U gives rise to a factor of 
^300 drop in the squared amplitude of this state after 
tunneling through the barrier. To trace the physical ori- 
gin of this substantial drop of transmission caused by 
the correlation effects in the FeOo.5 interface layer, we 
study the density of electronic states (DOS) of the sys- 
tem. As more than a half or this scattering state (54%) 
is contributed by Fe d z 2 orbitals, we concentrate on the 
behavior of Fe d z 2 states near the Fermi energy. Figure 
[3] shows the calculated DOS of d z 2 states of Fe atoms 
in the FeOo.5 interface layer. For the sake of compari- 
son we plot the calculated DOS curves for U = 4.6 eV 
and U = 0.0 eV (standard LSDA). We see a substantial 
difference in behavior of Fe d z 2 -electrons for "up" and 
"down" channels upon "turning on" the Coulomb repul- 
sion. For the spin-down states, the DOS at the Fermi 
level is practically unchanged, while for spin-up states it 
drops almost to zero when U = 4.6 eV is introduced. 
This is a well-known effect of LDA+ U : for the systems 
where correlation effects are strongly underestimated by 
standard LDA, the LDA+ U tends to split the DOS at 
the Fermi level into fully occupied and fully unnocupied 
peaks, significantly depleting DOS (or even forming the 
gap) at the Fermi levels. 

The performed DOS analysis explains well the trans- 
mission data, summarized in Table HI Taking into account 
the correlation effects in the FeOo.5 layer gives rise to a 
dramatic reduction of the electronic density of states at 
the Fermi level in the spin-up channel, which effectively 
"turns off" these electrons from the transmission process. 
In Table Q] this is reflected as a large drop in transmis- 
sion in the spin-up channel when a non-zero U is used 
in calculations. We decided to visualize this "turn-off" 




Energy, eV 

FIG. 3: Density of electronic states for the Fe electrons of 
d z 2 symmetry, calculated for the interface FeOo.5 layer. Both 
results for calculations with U = 0.0 eV and with U = 4.35 eV 
are presented. 
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FIG. 4: Real-space distribution of electronic density at the 
Fermi energy, calculated with U = 0.0 eV (left panel), and 
with U = 4.35 eV (right panel) for oxidated Fe-MgO-Fe sys- 
tem. Only states, having d z i symmetry, are taken into ac- 
count for Fe layers when calculating the density. 

culating their real-space distribution, and checking how 
this distribution is changed when the Coulomb repulsion 
is considered. In Figure 2] we plot the calculated real- 
space density of spin-up Fermi electrons for the interface 
part of the Fe/MgO/Fe system. We show only the elec- 
trons with d z 2 -symmetry of the wavefunction for all Fe 
atoms, while states with all symmetries are taken into ac- 
count for O and Mg atoms. This symmetry-based plot- 
ting allows us to trace the behavior only of those elec- 
trons which participate in the tunneling process. The 
picture shows that in the case of standard LSDA calcu- 
lations (left panel), the interface FeOo.s layer carries a 
significant density of d z i Fermi electrons, and it is even 
comparable to the electronic density in the deeper lay- 
ers of the Fe electrode. Moreover, these electrons couple 
with the oxygen p-electrons of the first MgO layer, cre- 
ating a non-zero density of Fermi electrons on this layer. 
When the on-site Coulomb repulsion is turned on, the 
d z 2 electrons with the Fermi energies are almost absent 
in the FeOo.s layer, as well as in the first layer of the 
MgO barrier. This picture illustrates that the proper 
consideration of the correlation effects in the interface 
FeOrj.5 layer leads to the localization of the Fe electrons 
with d z 2 symmetry, which does not allow them to partici- 
pate in the tunneling process and, as a result, is reflected 
as a dramatic reduction of the TMR characteristics of 
Fe/MgO/Fe system. 

D. Low-decay states and the symmetry of FeOo.s 
layer 

Let us now go back to the question of spectral com- 
position of the two low-decay scattering states, shown in 
Figl^l As we mentioned above, our calculations for the 



system with FeOo.s interface layer showed the presence of 
Fe d x 2_ y 2 orbitals in the low-decay scattering states: they 
contribute 30% to the first lowest-decay state (Figf2^) 
and are mostly forming (56%) the second one (Fig(2b). 
The orbitals of this symmetry do not contribute to the 
Ax-state in case of non-oxidated Fe/MgO/Fe interface 5 , 
and in the following we shall try to understand the phys- 
ical mechanisms of the influence of FeOo.s interface layer 
on the formation of the low-decay scattering states in 
Fe/MgO/Fe system. 

Butler et al* brought on the idea of "symmetry filter- 
ing" which appeared to be the central point of the ele- 
gant physical explanation of TMR effect in Fe/MgO/Fe 
trilayers. The authors demonstrated that Fe-bcc Bloch 
states of different symmetry show significantly different 
decay rates in the MgO barrier. The one with the small- 
est decay, which is able to couple to its counterpart in 
the opposite Fe lead, is that of Ai symmetry. This state 
transforms as a linear combination of s, p z , and d z 2 or- 
bitals. The analysis of Butler et al is based on a perfect 
cubic symmetry of Fe-bcc, matching also perfect cubic 
unit cell of MgO (although rotated by 45°). In the case 
of 100%-oxidated buffer layer (FeO) the situation does 
not change, as the point symmetry of all Fe sites is the 
same as in non-oxidated interface. However, the situation 
may change significantly if a low-symmetry FeOi-a; layer 
is formed at the Fe/MgO/Fe interface. In this case the 
admixing of other electronic states to Ai band may take 
place (strictly speaking, this low-decay band can not be 
classified as Ai anymore.). This picture correlates well 
with our results for the interface with FeOo.s iron oxide 
layer. Although the distribution of oxygen in this layer is 
still ordered, the symmetry of Fe sites is changed. There- 
fore our structure can be a good model to test the effects 
of symmetry changes in FeOi-n; layer on the composition 
of the low-decay scattering state. 

To clarify the role of symmetry of FeOo.s iron oxide 
layer in formation of the low-decay states in Fe/MgO/Fe 
interface, we now take a different approach and per- 
form a group theory analysis of the interaction of Fe 
lead states with the FeOi-a; layer. In case of 100%- 
populated FeO layer, all interface Fe atoms are occu- 
pying the positions with C<± v point symmetry group. In 
this group the <i-orbitals of iron and p-orbitals of oxygen 
of the first MgO layer transform according to the follow- 
ing irreducible representations: A\ => {d z 2,p z }; E => 
{ (d xz j dy Z 

),Px,P y }; B 1 =>■ {d xy }; B 2 =4> {d x 2^ y2 }. In 
this case the A\ orbitals are forming the Ai band. This 
band is the only one having the admixture of O p 2 -states, 
and has the lowest decay across the MgO barrier. The 
orbitals belonging to the two-dimensional representation 
E, form the A5 band. This band is doubly-degenerate 
and has a significantly higher decay in the MgO barrier. 
Finally, the orbitals belonging to B\ and B2 representa- 
tions form A2' and A2 bands, which are forbidden to in- 
teract with oxygen p-orbitals and show the highest decay 
in the MgO barrier. The results of this symmetry-based 
analysis are in excellent agreement with original physical 
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U=O.OeV U=4.6eV 




FIG. 5: The same quantity as in Fig(4] but calculated with 
contribution from Fe d x 2_ y 2 states only. 



picture presented by Butler et a&. 

Let us now perform the same analysis for the case of 
our FeOo.5 interface layer, where the point symmetry of 
Fe atoms is reduced. In this system the interface iron 
atoms have C 2v point symmetry, and the interacting Fe 
d- and O p-orbitals belong to the following irreducible 
representations^: A\ =>■ {d z 2 ,p z , d xy }; B\ {{d xz + 
d yz ),Px}\ B 2 {(d xz -d yz ),py}; A 2 =S> {d x 2^ y2 }. We 
immediately see the striking difference: the Ai represen- 
tation now includes both d z 2 and d xy orbitals while in 
the case of fully-symmetric FeO layer, the d xy orbitals 
were not allowed to interact with MgO p z -states by sym- 
metry rules. It is important to note here that due to the 
45° rotation of the Fe electrode's unit cell in our calcula- 
tion (as compared to the Fe-bcc unit cell) , the calculated 
d x 2_ y 2 orbitals found to contribute significantly to the 
low-decay scattering states, are in fact the d xy orbitals 
in Fe-bcc unit cell22. We see that our results of theoret- 
ical symmetry-based analysis excellently agree with the 
results of our numerical simulations. If we now speak 
the language of the symmetry bands, coming out of Fe- 
bcc (100) surface, we observe that introducing of ordered 
FeOo.5 interface layer in Fe/MgO/Fe interface leads to 
the mixing of Ai and A 2 r bands and their interaction 
with p 2 -states of oxygen in MgO layer. In this way, both 
Ai and A 2 < bands contribute now to the formation of 
low-decay scattering state in MgO barrier, which change 
significantly the TMR properties of the system. 

To complete our study, it is interesting to check if 
"turning on" the Coulomb repulsion influences interface 
Fe d x 2_ y 2 states in the same way as it collapses the d z 2- 
states studied above. We calculated the real-space distri- 
bution of electrons with the Fermi energy for both cases 
of U — and U = 4.35eV, where only the contribu- 
tion from d x 2_ y 2 orbitals is taken into account. The 
results are presented in FigJSJ We see that a non-zero 
Coulomb repulsion gives rise to a significant collapse of 



the d x 2_ y 2 states in the FeOo.5 layer. As a consequence, 
their coupling with oxygen states is dramatically reduced 
(notice the absence of charge in the first MgO layer for 
U = A.35eV case, Fig[5j right panel), which leads to the 
rf 3 .2_ 1 ,2-states being practically "switched off" from the 
tunneling process. It is important to notice that the same 
effect should also take place in non-oxidated system, or 
in the system with 100%-oxidated FeO layer, but in those 
cases the d x 2_ y 2 states are forbidden to interact with oxy- 
gen by symmetry and therefore this collapse should not 
influence in any way the tunneling characteristics of the 
interface. 

IV. CONCLUSION 

In conclusion, we presented the results of ab initio cal- 
culations of transport properties of Fe/MgO/Fe interface 
with FeOo.5 buffer layer. By combining NEGF-DFT for- 
malism with the LDA+ U approach, we explicitly took 
into account the on-site Coulomb repulsion in the Fe-0 
interface layer, where the [/-value was calculated from 
the first principles for a particular interface structure. 
Our results demonstrate a dramatic drop of TMR charac- 
teristic of the interface (down to ~60%), when Coulomb 
repulsion is taken into account in the FeOo.5 buffer layer. 
This reduction is found to be mostly due to a dramatic 
drop of transmission in spin-up channel for parallel con- 
figuration of iron leads. By analyzing the distribution of 
scattering states across the interface, as well as a real- 
space distribution of the Fermi electrons with a particu- 
lar symmetry of the wavefunction, we traced the physical 
origin of this transmission reduction. The results showed 
that for 50%-oxidated buffer layer (FeO . 5 ), this TMR 
drop is caused by a spacial collapse of Fe d z 2 and d x 2_ y 2 
electronic states at the Fermi level when non-zero [/-value 
is taken into account. Finally, we used the methods of 
group theory and studied the general type of interaction 
between MgO layer and the buffer FeOi-a, layer with dif- 
ferent point group symmetry. The results show that the 
reduction of the symmetry of the buffer FeOi-^ layer 
leads to the mixing of different symmetry bands of Fe- 
bcc bulk electrode, giving them the possibility to interact 
with MgO oxygen states and therefore to take part in the 
tunneling process. 

Appendix A: LDA+ U approach in MatDCal code 

The main idea of the LDA+ U method is to replace 
the LDA e — e Coulomb interaction energy of localized 
electrons in the energy functional of the system by a 
Hubbard-like termiiiSiSS: 

E LDA+U [p a (v),n a ] = E LDA [p a (r)} + Eu[n a ] - -Edcbv] 

( A1 ) 

Here Pa-(j) is a real-space electronic density for spin 
state a, Eld a \pa ( r )] is a standard LDA functional, and 
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n a is a density matrix of the subsystem of correlated 
orbitals. The Hubbard- type term of this modified func- 
tional can be written in the following way: 

Eu[n a ] =hj ^ »WW_o.+ 

rn,m>,„ (A2) 
+ ^{U ~ J) ^2 ri mG n m , 

Here we assume that the on-site Coulomb (U) and ex- 
change (J) interaction energies are independent of quan- 
tum numbers to. This expression also takes into account 
exchange: the interaction energy is U — J for electrons 
with the same spin a, while it is U for electrons with dif- 
ferent spin. The last term in (|A1|) is a so-called "double 
counting" correction whose role is to subtract from the 
functional the averaged electron-electron interaction in 
the correlated subsystem, which has already been taken 
into account within LDA: 



E dc [n„] = —N(N - 1) - 2 E N «( N ° - !) ( A3 ) 

a 

where N is the total number of electrons in correlated 
subsystem, and N a is the total number of electrons with 
particular spin a: 

N a =Y, n ™y, N = N a + N- a (A4) 

m 

Now, by inserting <fA2|) and ((A3]) in (|Al]) . and em- 
ploying the following expression for the total number of 
particles: 

m,cr m^m',m',cr m,m r ,er 

(A5) 



we arrive to a simple and computationally efficient ex- 
pression for LDA+ U functional: 

E LDA+U = E LDA + hp - J){N - O (A6) 

rn,(7 



Taking the functional derivative leads to the following 
expression for the orbital-dependent potential: 

V n L ° A+U = V LDA + hp J)(l - 2n ma ) (A7) 



As we notice, this form of LDA+ U potential de- 
pends on a single parameter which we call the "effective" 
Coulomb repulsion U e f / = U — J. This quantity can 
be calculated from the first principles for a given atomic 
system using a constrained LDA approach of Anisimov 
and Gunnarsso n 14 i 21 . 

To employ this scheme in localized basis set which the 
MatDCal package uses, we write the orbital-dependent 
part of the potential in operator form: 

V ma = V ma \ma)(ma\ (A8) 
where V ma is the second term in (|A7|) . and \ma)(ma\ is a 
projection operator on the atomic orbital basis function 
corresponding to the localized state, for which LDA+iJ 
treatment is applied. In this form we add the LDA+C/ 
potential in a straightforward manner to the MatD- 
CAL software where the Hamiltonian matrix elements 
are calculated. Having computed the matrix elements, 
we follow the standard computational scheme of NEGF 
method in localized basis set representation. The techni- 
cal details of this implementation can be found in several 
previous publication s 13 ' 16 ' 23 . 



* Electronic address: vladimir@physics.mcgill.ca 

1 S. Yuasa, T. Nagahama, A. Fukushima, Y. Suzuki, and 
K. Ando, Nat. Mater. 3, 868 (2004). 

2 S. Parkin, C. Kaiser, A. Panchula, P. Rice, B. Hughes, 
M. Samant, and S.-H. Yang, Nat. Mater. 3, 862 (2004). 

3 J. Moodera, L. Kinder, T. Wong, and R. Meservey, Phys. 
Rev. Lett. 74, 3273 (1995). 

4 A. Ney, C. Pampuch, R. Koch, and K. Ploog, Nature (Lon- 
don) 425, 485 (2003). 

5 W. Butler, X.-G. Zhang, and T. Schulthess, Phys. Rev. B 
63, 054416 (2001). 

6 D. Waldron, V. Timoshevskii, Y. Hu, K. Xia, and H. Guo, 
Phys. Rev. Lett. 97, 226802 (2006). 

7 J. C. R. P. G. Mather and R. A. Buhrman, Phys. Rev. B 
73, 205412 (2006). 

8 K. X. Y.Q. Ke and H. Guo, Phys. Rev. Lett. 105, 236801 
(2010). 



9 H. M. et al., Phys. Rev. Lett. 87, 076102 (2001). 

X.-G. Zhang, W. Butler, and A. Bandyopadhyay, Phys. 
Rev. B 68, 092402 (2003). 

1 G. M. et al, Phys. Rev. Lett. 100, 246803 (2008). 

2 S. Mirhosseini, K. Saha, A. Ernst, and J. Henk, Phys. Rev. 
B 78, 012404 (2008). 

3 J. Taylor, H. Guo, and J. Wang, Phys. Rev. B 63, 121104 
(2001). 

4 V. Anisimov and O. Gunnarsson, Phys. Rev. B 43, 7570 
(1991). 

5 P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka and J. 
Luitz, WIEN2k, An Augmented Plane Wave + Local Or- 
bitals Program for Calculating Crystal Properties (Karl- 
hcinz Schwarz, Techn. Universitt Wien, Austria), 2001. 
ISBN 3-9501031-1-2. 

6 Matdcal reference ... 

7 N. Troullier and J.L. Martins, Phys. Rev. B 43, 1993 



8 



(1991); K. Kleinman and D.M. Bylander, Phys. Rev. Lett. 
48, 1425 (1982). 

18 D. Waldron, P. Haney, B. Larade, A. MacDonald, and 
H. Guo, Phys. Rev. Lett. 96, 166804 (2006). 

19 V. Anisimov, J. Zaanen, and O. Andersen, Phys. Rev. B 
44, 943 (1991). 

20 V. Anisimov, I. Solovyev, M. Korotin, M. Czyzyk, and 
G. Sawatzky, Phys. Rev. B 48, 16929 (1993). 

21 G. Madsen and P. Novak, Europhys. Lett. 69, 777 (2005). 

22 Due to the semiinfinite nature of the iron lead, we assume 
the geometry of the electronic wavefunction, entering the 
FeOo.5 layer, to be the same as in the lead (Fe-bcc cu- 



bic lattice). However, the square lattice of FeOo.s layer, as 
well as the MgO lattice, are 45° rotated with respect to 
Fe-bcc unit cell. This fact is taken into consideration when 
we classify Fe states, coming out of the lead, in the new 
symmetry of FeOo.5 and MgO layers. The most important 
in this context is that the d x i_ y i orbital, which significant 
contribution to the scattering states was obtained by nu- 
mericall computations, is in fact the d xy orbital in Fe-bcc 
lattice. 

E. Marcotte, M.Sc. thesis, McGill University. (2009). 



